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AIRBLAST  SIMULATIONS  USING  FLUX-CORRECTED  TRANSPORT  CODES 


INTRODUCTION 

The  method  of  Flux-Corrected  Transport  (FCT)  was  developed  about  ten 
years  ago  (Boris  and  Book,  1973;  Book,  et  al,  1975;  Boris  and  Book,  1976)  as 
a  means  for  solving  systems  of  hyperbolic  equations  in  such  a  way  that 
physically  positive  quantities  remain  positive.  Its  principal  application 
has  been  to  compressible  fluids,  i.e.,  fluids  in  which  some  of  the  flow 
speeds  are  comparable  with  or  greater  than  the  local  speed  of  propagation  of 
waves  in  the  medium.  Examples  include  plasmas,  combusting  systems,  and  gas- 
dynamic  flows. 

Subsequently  the  method  has  been  extended  and  generalized  in  a  number  of 

important  respects.  These  include  the  construction  of  general-purpose 

algorithms  for  solving  fluid  equations  on  moving  nonuniform  grids  in  a 

curvilinear  coordinate  system  (Boris  1976),  strictly  multidimensional  FCT 

algorithms  (Zalesak  1979)  (as  opposed  to  timestep-splitting  of  one- 

dimensional  algorithms),  and  variants  on  these.  In  addition,  numerous 

techniques  have  been  developed  by  other  workers  which  incorporate  the 

positivity-preserving  property  in  other  ways  (Harten,  1978;  Van  Leer  1979; 

Collela,  1983).  The  weight  of  opinion  among  computational  scientists  who 

deal  with  problems  of  compressible  flow  now  strongly  affirms  the  superiority 

of  positivity-preserving  methods  over  conventional  ones.  In  a  few  words, 

advantages  are  improved  robustness,  flexibility,  resolution,  and  freedom  from 

nonphysical  numerical  artifacts. 

Manuacript  approved  February  2,  1984. 


In  this  paper  we  describe  some  applications  of  FCT  to  airblast  problems, 
i.e.,  to  calculations  of  the  blast  wave  and  fireball  produced  by  an  explosion 
in  the  atmosphere  near  the  ground.  As  this  is  one  of  the  applications  for 
which  FCT  was  originally  devised,  it  is  not  surprising  that  the  method  yields 
satisfactory  results  in  a  wide  range  of  cases  with  very  little  need  for 
computational  development.  Nevertheless,  some  of  the  refinements  we  have 
employed  in  this  connection  are  of  interest  both  in  their  own  right  and  for 
other  users  who  may  need  to  carry  out  similar  calculations. 

In  most  of  the  airblast  calculations  we  have  carried  out,  four  distinct 
stages  are  discernible:  (i)  initialization,  usually  employing  a  Chapman- 
Jouguet  model  (Kuhl  et  al,  1981)  to  describe  conditions  at  the  time  the 
detonation  wave  reaches  the  periphery  of  the  explosive  charge;  (ii)  free- 
field  (one-dimensional)  evolution,  before  the  shock  reaches  the  ground;  (iii) 
subsequent  two-dimensional  evolution  using  axisymmetric  (r-z)  geometry;  and 
(iv)  postprocessing  of  the  resulting  solution  to  analyze  and  display  the 
results. 

The  advances  described  in  this  paper  relate  to  all  but  the  first  of 
these  states.  They  may  be  cataloged  as  follows:  refinements  in  the  free- 
field  solution;  development  of  a  regridding  algorithm  for  the  early-time 
(high-overpressure)  two-dimensional  solution;  improvements  in  boundary 
conditions  in  the  two-dimensional  stage;  and  exploitation  of  tracer 
particles  as  a  diagnostic,  particularly  in  the  late  stages  of  fireball 
evolution,  after  shock  breakaway.  In  the  body  of  the  paper  one  section  is 
devoted  to  each  of  these  topics. 
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REFINEMENTS  IN  THE  FREE-FIELD  SOLUTION 

Figure  1  shows  one-dimensional  (spherical)  profiles  calculated  using  a 
form  of  the  general-purpose  FCT  algorithm  ETBFCT  (Boris  (1976).  Although 
these  profiles  are  accurate  and  physically  correct  on  the  average,  they 
exhibit  numerous  flat  shelves  and  accompanying  sharp  drops,  especially  in  the 
region  of  the  rarefaction  wave.  This  phenomenon  is  known  as  "terracing.*'  In 
this  instance  we  attribute  its  presence  to  the  nonlinear  action  of  the  flux 
limiter  stage  of  FCT  on  dispersive  errors  from  the  transport  stage. 

To  explain  what  this  means,  let  us  describe  a  class  of  FCT  algorithms 
which  includes  ETBFCT  as  a  particular  case.  To  advance  the  one-dimensional 
continuity  equation 


3p  3  . 

It  *  “  (PV) 


(1) 


one  time  step  on  a  uniform  mesh  with  a  given  flow  velocity  v,  we  suppose  the 

old  cell-centered  values  p°  and  v^  are  known.  A  general  three-point 

o 

conservative  finite-difference  operator  acting  an  p^  can  be  written  in  the 
form 

To 

Pj  "  Pj  '  £j+l/2  Pj+L/2  +  Ej-1/ 2  Pj-l/2 


(2) 


+  Yj+l/2  (pj+l  "  Pj}  Vl/2  (Pj  "  PJ-1)* 


where  Pj+j/2  *  ~2  Pj+  Pj+1^*  To  °^tain  a  first-order  approximation  to  Eq. 

(I),  we  must  define  ej+1/2  to  t*ie  Courant  number  vj+i/2^t^X>  where 

vj+l/2  m\  ^vj  +  Vj+P’  likewise,  for  a  conventional  scheme  with  second- 

1  2 

order  (spatial)  accuracy,  we  must  let  y  ^ ^  “  *£  ej+i/2*  In  FCT»  however,  a 

T 

numerical  diffusion  is  applied  to  p^  according  to 
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(3) 


TD  T  ,  o  o.  .00. 

P<  ’  0-.  +  V1+l/2Cej+l  -  V  '  vj-l/2(pj  -  Vi>  • 

TD 

Evidently  the  y's  and  the  v's  combine  in  the  same  way  in  p^  ,  so  we  leave 
both  unspecified  for  the  moment. 

If  we  now  applied  another  diffusion  operation,  this  time  with  the 

opposite  sign,  it  would  just  tend  to  cancel  the  diffusion  terras  already 

TD  n 

present  in  .  Instead  we  write  for  the  new  densities  p^ 


n  TD  ,c  ,  c 
Pj  "  Pj  "  *3+1/2  +  * j — 1  /  2  * 


(4) 


where  the  "corrected"  fluxes  differ  from  the  "raw”  fluxes  A.,.,.  = 

J+l/2  j+1/2 

Uj+i/2  ^pj+i  ~  Pj)  in  two  ways:  we  use  {pT }  instead  of  {p°}  in  the  definition 
of  the  raw  fluxes  (this  allows  us  to  define  Y  so  as  to  optimize  some  property 
in  the  difference  scheme,  as  will  be  seen  shortly);  and  we  correct  the 

fluxes,  so  that  the  antidiffusion  process  (which  evidently  tends  to  make  all 

TD 

gradients  steeper)  can  enhance  no  extrema  already  present  in  {p^  },  nor 
introduce  any  new  ones.  The  simplest  formula  for  achieving  this  in  all 
possible  situations  is  that  used  in  "strong  flux  limiting:” 


<t>j+l/2  -  S  •  {max[0,  min  C  |  +j+1/2 1 »  Aj-i/2’  Aj+3/2)]^  •  (5) 

TD  TD 

where  S  -  sign  4>j+1/2  and  Aj+2/2  =  Pj+J  -  Pj  . 

At  most  points  in  a  gently  varying  profile  no  correction  is  required  and 
c 

^j+1/2  *  ^j+1/2*  *n  t*lat  case  we  can  perform  a  Von  Neumann  analysis, 
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calculating  the  complex  propagator  or  amplification  factor  A  **  Ar  +  LA^  - 
pj/pj  for  a  sinusoidal  density  profile  p°  *  exp(ijk  <5x),  where  k  is  the  wave 
number.  Writing  0  =  k6x  and  e  *  v5t/  &c,  we  expand  in  powers  of  0  the 
amplification 


|a|  *  1  +  a2  02  +  a4  04  +  ... 


and  the  relative  phase  error 


R  -  tan-1  <w 


-  1  *  R2  02  +  R4  04  +  ., 


It  is  easy  to  show  that  the  condition  that  A2  vanish  is 


1  o 

Y+v—  u-  —  ez. 


The  condition  that  R2  vanish  is 


(6) 


(7) 


(8) 


U 


(9) 


There  remains  one  free  parameter  which  can  be  chosen  so  as  to  build  some 
desirable  property  into  the  algorithm  (this  is  the  reason  for  introducing  y). 
The  choice  Y  “  e2  confers  no  particular  benefit.  Letting  y  5  0  as  in 

ETBFCT  (Boris,  1976)  reduces  the  operation  count  and  yields  a  useful  general- 

purpose  algorithm,  which  is,  however,  subject  to  the  terracing  problems 

mentioned  previously.  There  are  two  other  obvious  candidates:  y  “  4-  e2 

4 

which  implies  A^  ■  0,  and  y  m  e2,  which  implies  R^  -  0. 
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Figure  2  shows  this  for  a  test  involving  passive  advection  with  a 
uniform  velocity.  A  test  profile,  initially  consisting  (Fig.  2a)  of  a  half 
ellipse,  connected  by  two  quarter  ellipses  to  a  straight  line  of  zero  slope, 
is  advected  with  a  Courant  number  e  =  0.2.  The  results  obtained  using  the 
general-purpose  algorithm  (Fig.  2b)  and  the  new  multicoefficient  versi.n 
(Fig.  2c)  are  shown  for  comparison.  It  is  clear  that  the  latter  represents  a 
distinct  improvement,  a  conclusion  that  continues  to  hold  for  other  test 
profiles. 

The  results  of  a  different  test  are  shown  in  Fig.  3.  In  Fig.  3a,  the 
profile  obtained  from  an  analytic  solution  of  the  bursting  diaphragm 
(Riemann)  problem  with  10-to-l  density  ratio  and  103-to-l  pressure  are 
plotted,  while  in  Fig.  3b  those  calculated  with  the  new  FCT  algorithm  are 
shown  for  comparison.  Here  again  the  results  agree  fairly  closely.  The 
agreement  is  worst  near  the  contact  surface,  which,  unlike  the  shock,  has  no 
physical  self-steepening  mechanism. 

In  chemical  explosions  a  characteristic  feature  is  the  contact  surface 
between  the  detonation  products  and  the  air  compressed  ahead  of  the  blast 
(Fig.  4).  In  the  absence  of  diffusion  (usually  a  good  approximation  in  an 
actual  explosion,  if  not  numerically),  the  contact  surface  remains  sharp,  and 
the  two  media  differ  in  density  across  this  discontinuity  by  ~  30%. 

In  an  airblast,  a  layer  of  air  is  compressed  beneath  the  burst,  so  that 
the  detonation  products  never  reach  the  ground.  The  shock  wave  propagates 
down  to  the  ground,  reflects  upward,  and  propagates  back  through  the  layer  of 
compressed  air.  When  it  encounters  the  contact  surface,  it  is  partly 
transmitted  and  partly  reflected  downward  again.  This  process  can  repeat 
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several  times.  At  finite  ground  range  the  shocks  reflect  obliquely  and  soon 
propagate  away  from  the  fireball,  but  directly  beneath  the  burst  the  process 
is  only  terminated  when  the  rising  fireball  recedes  from  the  ground 
sufficiently  far.  (Experimentally  and  computationally  the  multiply  reflected 
shocks  usually  become  indistinguishable  from  noise  after  the  second  or  third 
reflection. ) 

Therefore  it  is  important  computationally  to  model  the  contact  surface 
accurately  if  these  shock  reverberations  are  to  be  described  properly. 

If  the  contact  surface  in  the  simulation  is  diffuse,  the  shock  reflected  from 
it  will  be  weak  or  undetectible.  (If  physical  diffusion  processes  were 
present  experimentally  but  absent  in  the  simulation,  however,  the  latter 
would  predict  excessively  strong  reflection.)  The  results  shown  in  Fig.  4a, 
calculated  using  200  zones,  while  close  to  the  physically  correct  profiles, 
exhibit  considerable  diffusion  at  the  contract  surface.  This  is  in  contrast 
with  the  results  shown  in  Fig.  4b,  obtained  using  400  zones,  in  which  the 
contact  surface  is  much  more  distinct.  A  consequence  of  improving  the 
resolution  is  shown  in  Fig.  5,  which  depicts  the  pressure  history  that  a 
sensor  beneath  the  bursts  shown  in  Fig.  4  would  measure.  The  peaks  are  due 
to  the  initial  shock  and  the  shock  re-reflected  from  the  contact  surface, 
respectively.  Note  that  both  maxima  increase  when  the  resolution  improves, 
but  the  second  peak  rises  more.  Although  in  this  case  further  enhancement  of 
the  diffusion  makes  little  difference,  it  is  clear  that  the  second  peak  might 
in  some  cases  be  raised  above  the  first.  Thus  a  qualitative  change  in  the 
solution  can  result  from  improving  the  resolution.  This  is  of  course  also 
applied  to  the  accuracy  of  the  two-dimensional  results  calculated  after  the 
one-dimensional  solution  is  laid  down  on  2D  mesh. 
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REGRIDDING  ALGORITHM  FOR  TWO-DIMENSIONAL  SOLUTION 


To  study  the  transition  to  Mach  reflection  extremely  high  resolution  is 
required,  as  may  be  seen  from  the  following  discussion.  The  Mach  stem  at 
transition  has  vanishing  length.  If  we  are  restricted  to  resolution  of 
features  several  (say,  a  factor  of  three)  times  the  mesh  size  5x  ~  6y,  then 
the  ground  range  at  which  transition  occurs  will  be  overestimated  by  an 
amount  Si  equal  to  this  length  divided  by  the  tangent  of  the  opening  angle  x> 

which  is  typically  a  small  fraction  of  a  radian,  so  the  6>Z  ~  30  6x.  Since 

core  limitations  restrict  us  to  meshes  with  a  total  of  no  more  that  104  -  105 
zones  and  since  transition  typically  takes  place  at  a  ground  range  comparable 

with  the  height  of  burst  (HOB),  if  we  use  a  uniform  mesh  the  error  <$£  would 

be  as  much  as  ten  or  twenty  percent  of  the  HOB. 

Fortunately,  this  is  unncessary.  We  can  grid  up  near  where  transition 
takes  place,  using  cells  smaller  by  one  or  two  orders  of  magnitude  in  each 
direction,  with  a  proportionate  reduction  in  61.  Since  the  mesh  is 
rectilinear,  this  implies  the  presence  of  extremely  elongated  (non-square) 
zones  on  portions  of  the  grid,  with  consequent  errors  in  the  solution.  This 
creates  no  difficulties  if  the  details  of  the  flow  there  are  unimportant  and 
if  those  regions  do  not  generate  errors  which  propagate  into  the  regions  of 
chief  interest. 

However,  a  shock  wave  propagating  from  coarsely  gridded  region  abruptly 
into  a  finely  gridded  one  sheds  sound  waves  as  it  steepens  up,  seriously 
degrading  the  solution  at  and  behind  the  front  (Boris  and  Book  (1976)). 

Hence  it  is  necessary  to  ensure  that  the  shock  at  the  front  of  the  blast  wave 
be  allowed  to  propagate  some  distance  through  a  finely  gridded  region  prior 
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to  reaching  the  point  of  reflection,  that  the  transition  be  gradual,  and  most 
significantly,  that  this  be  true  not  just  at  transition  but  at  all  times  up 
until  then.  The  only  way  to  avoid  using  an  inordinate  number  of  fine  zones 
is  to  move  the  highly  resolved  region  with  the  reflection  point,  starting  at 
the  instant  the  blast  wave  reaches  the  ground  and  continuing  as  long  as  the 
Mach  reflection  process  is  the  interest  (Fig.  6).  This  method  exploits  the 
"sliding  rezone"  (continuous  regridding)  capability  of  the  ETBFCT  algorithm. 

Typically  we  start  with  a  mesh  ~  150  x  150  zones,  having  ~  100  fine 
zones  of  dimension  6x  =  Sy  and  ~  50  coarse  zones  of  dimension  Ax  =  Ay  in  each 
coordinate  direction,  with  Ax/  <5x  =  10.  The  fine  zones  are  those  closest  to 
the  axes.  Let  the  cell  boundaries  be  specified  by  {x^  }  and  {y ^  }.  We 
introduce  a  smooth  transition  region  by  diffusing  these  interfaces  repeatedly 
according  to 


xi  '  \  +  "<x1+i 


*1  +  Vi> 


(10) 


'  ’  +  n(Vi 


2yj  + 


(ID 


where  the  diffusion  coefficient  n  is  of  order  0.1.  The  diffusion  process  is 
repeated  N  times  (N  ~  30)  to  generate  a  transition  region  N  cells  wide.  The 
initial  reflection  necessarily  occurs  in  the  fine-zoned  region. 

As  the  reflection  point  moves  outward,  the  grid  is  redefined  each  cycle. 
We  search  on  density  (or  pressure)  to  find  the  smallest  radial  displacement 
such  that  conditions  deviate  from  ambient  by  a  significant  amount  (3%). 
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Using  this  as  a  center,  we  define  new  cell  boundaries,  smooth  {x^  }  and 
{y^  }  as  before,  and  deposit  the  transported  fluid  quantities  into  the  cells 
of  this  new  grid  at  the  conclusion  of  the  transport  calculation.  As  the 
fine-zoned  region  marches  outward,  coarse  cells  disappear  ahead  of  it  and 
reappear  at  small  radii  (Fig.  6). 

This  method  has  proved  successful  on  a  number  of  HOB  calculations.  A 
certain  amount  of  tuning  is  necessary  to  ensure  that  enough  zones  are  present 
in  both  the  highly  resolved  and  transition  regions,  but  there  is  no 
difficulty  in  tracking  the  reflection  point,  nor  is  any  significant  detail 
lost  through  diffusion. 

BOUNDARY  CONDITIONS 

The  conditions  for  reflection  at  a  planar  solid  boundary  in  a  perfectly 
invisid  flow  are 


3p  9p  3vt 

3n  *  dn  3n  Vn 


(12) 


where  represents  a  normal  derivative  and  the  subscripts  n  and  t  label  the 

normal  and  tangential  components  of  velocity.  These  conditions  are  imposed 

at  the  interface  between  the  last  cell  of  the  computational  mesh  and  a  guard 

cell  juxtaposed  to  it.  To  ensure  that  the  normal  derivative  of  a  quantity  f 

vanish,  it  is  sufficient  to  set  the  value  in  the  guard  cell  equal  to  that  in 

the  last  internal  cell,  i.e.,  fN+^  =  f^.  To  ensure  that  the  quantity  itself 

=  -f„.  (Special  attention  must  be  paid  to  the 
N 


vanish,  we  simply  set  f 


equation  describing  the  component  of  momentum  normal  to  the  boundary  in  order 


that  no  momentum  be  diffused  through  the  interface.)  Equation  (12)  applies 
to  any  symmetry  plane  or  axis,  and  so  to  the  axis  of  symmetry  in  the  HOB 
problem  at  r=0.  At  large  values  of  r  and  z  (i.e.,  at  the  peripheral  boundary 
—  the  outside  and  top  of  the  mesh),  we  impose  the  ambient  conditions  p  ■ 

Pa.  P  -  Pa,  v  -  0. 

This  is  satisfactory  until  shocks  propagate  to  the  boundary.  If  we  are 
interested  in  following  the  shocks  further,  it  is  necessary  at  this  point  to 
begin  moving  the  grid  outward  (as  described  in  the  previous  section),  or  to 
lay  down  the  solution  on  a  larger,  coarser  mesh.  If  we  are  interested  in 
characteristic  late-time  phenomena  such  as  vortex  formation  and  fireball 
rise,  however,  shocks  are  irrelevant  and  it  is  appropriate  to  use  simple 
outflow  (or  transmission)  boundary  conditions,  namely 


3p  3p  3v 
3n  =  3n  3n 


(13) 


Higher-order  extrapolation  schemes  have  been  recommended  in  the 
literature  (see  e.g.  ,  Roache  1972),  but  have  not  been  found  particularly 
useful  in  connection  with  FCT,  apparently  because  of  the  way  they  interact 
with  FCT's  tendency  to  clip  off  the  tops  of  peaked  profiles.  At  very  late 
times  a  small  but  significant  modification,  suggested  by  Boris  (personal 
communication,  1983)  has  been  employed.  This  involves  introducing  a 
relaxation  term  at  the  peripheral  boundary  in  the  form 


’  N+l 


f N  +  °(fa  "  V 


(14) 
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Here  the  relaxation  coefficient  satisfies  a  *  c  dt/L,  where  c  is 
the  ambient  speed  of  sound  and  L  is  a  typical  scale  size  (e.g.,  fireball 
cloud  diameter)  in  the  problem.  Typically  a  is  small,  a  ~  10“  2  -  ID-3.  Its 
effect  is  to  cause  the  peripheral  values  of  the  fluid  variables  to  relax 
toward  their  ambient  values,  simulating  the  influence  of  inflow  from  the 
surrounding  atmosphere. 

Also  at  very  late  times  (typically  tens  of  seconds),  gravity  must  be 
included  in  the  vertical  momentum  equation  to  describe  the  effect  of 
buoyancy.  The  equation  we  solve  has  the  form 


3(  Pv) 
3t 


+  <r(Xlv)  +  (pv2)  (p  " 


dz 


dz 


Pa) 


(15) 


+  ( p  -  p  )g  «  0. 

a 

The  ambient  quantities,  which  are  taken  from  tables,  are  the  same  as 
those  used  in  the  initial  conditions.  These  represent  a  hydrostatic  balance 
which  truncation  errors  will  not  destabilize. 

The  energy  equation  also  acquires  an  additional  term  due  to  gravity; 
this,  however,  plays  only  a  minor  role. 

Figure  7  shows  the  late  time  evolution  of  a  fireball  calculated  for  a 
600-ton  burst  of  high  explosive  (HE)  set  off  at  a  height  of  50.6  m.  Note  the 
elongation  of  the  cloud  and  the  characteristic  vortices  which  appear 
especially  clearly  in  the  HE  density  contours. 


GRAPHICS  DEVELOPMENTS 


Approximately  one  quarter  of  our  computational  costs  are  accumulated  in 
generating  pictures  showing  the  results  of  the  calculation.  These  take  the 


form  of  pressure  histories;  plots  of  the  surface  pressure  and  density  vs 
range;  velocity  vector  (arrow)  plots;  contours  of  constant  total  and  HE 
product  mass  density,  total  and  internal  energy  density,  pressure,  and 
vorticity;  instantaneous  locations  or  trajectories  of  passively  advected 
tracer  particles;  and  movies  of  the  time  variation  of  these  displays.  Both 
the  latter  and  the  use  of  distinguishing  colors  are  particularly  revealing, 
although  examples  can  not  be  presented  here.  In  three-dimensional  work  we 
have  carried  out  (Fry  and  Book  1983,  unpublished),  the  problems  of  displaying 
results  become  significantly  greater. 

The  approach  we  have  taken  is  to  perform  the  calculation  in  its  entirety 
first,  often  with  several  restarts,  producing  only  enough  pictures  to  make 
sure  everything  is  going  well,  and  dumping  out  the  results  (essentially  the 
two-dimensional  p,  u,  v,  and  p  arrays)  at  intervals  onto  magnetic  tape.  This 
means  that  If  we  require  a  posteriori  a  quantity  which  cannot  be  deduced  from 
the  data  in  the  dump  files  (such  as  the  total  impulse  on  the  ground  sensors 
or  the  pressure  history  at  a  location  where  no  sensor  existed),  the 
calculation  has  to  be  rerun.  This  is  true  also  if  the  dump  interval  chosen 
is  too  long,  but  In  many  cases  intermediate  values  can  be  obtained 
satisfactorily  by  interpolation.  For  example,  a  movie  of  some  700  frames  was 
made  using  linear  interpolation  between  twenty-seven  dumps.  (Cubic 
interpolation  would  have  required  generating  a  new  data  type  containing  only 
the  arrays  being  processed,  since  the  contents  of  four  dumps  cannot  be  fitted 
into  core  simultaneously.)  The  graphics  calculations,  like  the  hydrodynamics 
calculation  themselves,  were  performed  on  a  Texas  Instruments  ASC,  a  two-pipe 
vector  machine  with  approximately  half  a  million  words  of  core.  The  plots 
were  output  on  a  Tektronix  4014  scope;  the  color  movie  was  generated  using  a 
dicomed  graphics  system.  This  approach,  in  which  we  generated  graphics  by 
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postprocessing  a  limited  number  of  dump  files,  is  a  compromise  between 
generating  all  the  graphics  in  the  course  of  the  calculation  and  saving  the 
results  from  every  timestep.  (This  would  entail  prohibitive  disk  or  I/O 
requirements.)  In  general  it  has  worked  well  particularly  as  we  have  learned 
from  experience  what  quantities  should  be  saved  and  how  often  they  should  be 
dumped . 

In  concluding,  we  remark  that  passive  advection  of  tracer  particles  is 

among  the  most  useful  of  the  graphical  display  we  have  used,  particularly  in 

movies.  Figure  8  shows  tracer  particle  histories  over  two  time  intervals  for 

the  calculation  described  previously.  Note  that  four  distinct  vortices  are 

revealed  in  these  pictures,  two  in  the  "forward”  direction  and  two  reversed. 

In  previous  work  (both  experimental  and  in  simulations),  only  one  of  each  had 

been  detected.  Moreover,  we  are  not  restricted  to  passive  advection,  which 

is  appropriate  to  the  motion  of  infinitesimal  particles  in  the  blast  wind 

fields.  We  have  also  modeled  finite-sized  particles  (dust  and  debris)  by 

balancing  inertial  forces  against  gravity  and  the  Stokes  drag  according  to 

dV 

m-j^-  -  -mg  +  D(v  -  V),  (16) 

where  m  and  V  are  the  particle  mass  and  velocity,  respectively,  and  D  is  the 
drag  coefficient.  By  using  (16)  with  a  spectrum  of  particle  masses  {m^ } 
and  appropriate  initial  conditions  obtained  from  an  ejecta  model,  we  can 
generate  realistic  maps  showing  how  must  material  is  entrained  by  the  blast 
and  where  it  is  located  at  any  given  time.  This  is  invaluable  for  comparison 
with  observations  of  field  tests,  where  the  obscuring  debris  clouds  are  among 
the  most  conspicuous  features. 
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P.O.  Box  3707 

Seattle,  Washington  98124 

01CY  Attn  S  Strack 

01CY  Attn  Aerospace  Library 

01CY  Attn  M/S  42/37  R  Carlson 

California  Research  &  Tech  Inc. 

4049  First  Street 
Livermore,  CA  94550 
01CY  Attn  D  Orphal 


Applied  Theory,  Inc. 

1010  West  Wood  Blvd 

Los  Angeles,  CA  90024 

(2  Cys  if  unclass  or  1  Cy  if  Class) 

01CY  Attn  J  Trulio 
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Astron  Research  &  Engineering 
1901  Old  Middlef ield  Way  #15 
Mountain  View,  CA  94043 
01CY  Attn  J  Huntington 


APTEC  Associates  Inc. 

26046  Eden  Landing  Road 
Hayward,  CA  94545 
01CY  Attn  S  Gill 

AVCO  Research  &  Systems  Group 
201  Lowell  Street 
Wilmington,  MA  01887 
0 ICY  Attn  Library  A830 


ACUREX  Corp. 

485  Clyde  Avenue 
Mountain  View,  CA  94042 
01CY  Attn  C  Wolf 


Agbabian  Associates 
250  N  Nash  Street 
El  Segundo,  CA  90245 
01CY  Attn  M  Agbabian 

Applied  Research  Associates,  Inc. 
.2601  Wyoming  Blvd  NE  Suite  H-l 
Albuquerque,  NM  87112 
01CY  Attn  J  Bratton 
01CY  Attn  N  Higgins 

Department  of  the  Interior 
Bureau  of  Mines 

Bldg  20,  Denver  Federal  Center 
Denver,  CO  80225 
(Unclasified  Only) 

01CY  Attn  Tech  Library 


Sandia  National  Lab 

P.0.  Box  5800 

Albuquerque,  NM  87185 

All  Class  Attn  Sec  Control  0FC  FOR 

01CY  Attn  A  Chaban 

0 ICY  Attn  L  Hill 

01CY  Attn  Org  1250  W  Brown 

01CY  Attn  A  Chabia 

01CY  Attn  W  Roherty 

01CY  Attn  3141 

01CY  Attn  L  Vortman 

01CY  Attn  J  Banister 


BDM  Corp. 

7915  Jones  Branch  Drive 
McLean,  VA  22102 
01CY  Attn  A  Lavagnino 
01CY  Attn  T  Neighbors 
01CY  Attn  Corporate  Library 

Aerospace  Corp. 

P.0.  Box  92957 

Los  Angeles,  CA  90009 

01CY  Attn  H  Mirels 

OICY  Attn  Technical  Info  Services 

Analytic  Services,  Inc. 

400  Army-Navy  Drive 
Arlington,  VA  22202 
01CY  Attn  G  Hesselbacher 

Central  Intelligence  Agency 
Washington,  DC  20505 
01CY  Attn  0SWR/NED 


Director 

Federal  Emergency  Management  Agency 

Natinal  Sec  Ofc  Mitigation  &  Rsch 

1725  I  Street  NW 

Washington,  DC  20472 

(All  class  attn  B105  DOC  Control  FOR 

01CY  Attn  Mitigation  &  Rsch  Div. 

Los  Alamos  National  Scientific  Lab 

Mail  Station  5000 

P.0.  Box  1663 

Los  Alamos,  NM  87545 

Classified  only  to  mail  station  5000 

01CY  Attn  R  Whittaker 

01CY  Attn  C  Keller 

01CY  Attn  M.T.  Sanford 

01CY  Attn  MS  364  Class  Reports  LIB 

01CY  Attn  E  Jones 
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Lovelace  Biomedical  &  Environment 
Rsch  Institute,  Inc. 

P.0.  Box  5890 

Albuquerque,  NM  87115 

01CY  Attn  R  Jones  (unci  only) 


Sandia  Laboratories 
Livermore  Laboratory 
P.O.  Box  969 
Livermore,  CA  94550 


Department  of  Energy 
Washingotn,  Dc  20545 
01CY  Attn  OMA/RD&T 


Lawrence  Livermore  National  Lab 

P.O.  Box  808 

Livermore,  CA  94550 

01CY  Attn  L-90  R  Dong 

01CY  Attn  L-250  J  Hearst  (Class  L-203) 

01CY  Attn  L-90  D  Norris  (Class  L-504) 

01CY  Attn  L-7  J  Kahn 

01CY  Attn  D  Glenn 

01CY  Attn  L-437  R  Schock 

01CY  Attn  Technical  Info  Dept  Library 

01CY  Attn  L-200  T  Bltkovich 


Assistant  Chief  of  Staff 
Intelligence 

Department  of  the  Air  Force 
Washington,  DC  20330 
01CY  Attn  In  Rm  4A932 


Oak  Ridge  National  Laboratory 

Nuclear  Division 

X-10  Lab  Records  Division 

P.O.  Box  X 

Oak  Ridge,  TN  37830 

01CY  Attn  Civil  Def  Res  Proj 

01CY  Attn  Att  Central  Rsch  Library 

Department  of  Energy 
Albuquerque  Operations  Office 
P.O.  Box  5400 
Albuquerque,  NM  87115 
01CY  Attn  CTID 

Department  of  Energy 

Nevada  Operations  Office 

P.O.  Box  14100 

Las  Vegas,  NV  89114 

01CY  Attn  Mail  &  Records  for 

Technical  Library 

Commander 

Rome  Air  Development  Center  AFSC 
Griffiss  AFB,  NY  13441 
(Desires  No  CNWDI) 

01CY  Attn  TSLD 

Strategic  Air  Command 

Department  of  the  Air  Force 

Offutt  AFB,  NB  68113 

01CY  Attn  NRI-Stinf o  Library 

01CY  Attn  XPFS 

01CY  Attn  Int  J  McKinney 

Director 

Air  University  Library 
Department  of  the  Air  Force 
Maxwell  AFB,  AL  36 1 12 
(Desires  NO  CNWDI) 

01CY  Attn  Aul-LSE 

Assistant  Chief  of  Staff 
Studies  &  Analyses 
Department  of  the  Air  Force 
Washington,  DC  20330 
01CY  Attn  AF/SAMI  (Tech  Lib) 
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Assistant  Secretary  of  the  AF  Air  Force  Geophysics  Laboratory 

Research,  Development  and  Logistics  Hanscom  AFB,  MA  01731 

Department  of  the  Air  Force  01CY  Attn  LWW  K  Thompson 

Washington,  DC  20330 

01CY  Attn  SAFALR/Dep  for  Strat  &  Space  Sys 


Air  Force  Institute  of  Technology 
Air  University 

Wright-Patterson  AFB,  OH  45433 
(Does  Not  Desire  Classified  Documents) 
01CY  Attn  Library 

Air  Force  Weapons  Laboratory,  AFSC 

Kirtland  AFB,  NK  87117 

01CY  Attn  NTES-C  R  Henny 

01CY  Attn  NTED-I 

01CY  Attn  NTED  R  Matalucci 

01CY  Attn  NTF  M  Plamondon 

01CY  Attn  NT  D  Payton 

01CY  Attn  NTED-A 

0 ICY  Attn  NTES-G 

01CY  Attn  SUL 

01CY  Attn  DEX 

01CY  Attn  NTES-S 

01CY  ATtn  NTEO 

0  ICY  Attn  DEY 


Deputy  Chief  of  Staff 
Logistics  &  Engineering 
Department  of  the  Air  Force 
Washington,  DC  20330 
01CY  Attn  LEEF 

Office  of  the  Chief  of  Naval  Operations 
Washington,  DC  20350 
01CY  Attn  OP  981 
01CY  Attn  OP  03EG 

Commander 

Naval  Surface  Weapons  Center 
Dahlgren,  VA  22449 

01CY  Attn  Tech  Library  &  Info  Svcs  Br 
Commander 

Naval  Weapons  Center 
China  Lake,  CA  93555 
01CY  Attn  Code  3201  P  Cordle 
01CY  Attn  Code  266  C  Austin 
01CY  Code  233  (Tech  Lib) 


Headquarters 

Air  Force  Systems  Command 
Andrews  AFB,  DC  20334 
01CY  Attn  DLWM 
01CY  Attn  DLW 

Ballistic  Missile  Office/MN 
Air  Force  Systems  Command 
Norton  AFB,  CA  92409 
(Minuteman) 

01CY  Attn  MNNXH  G  Kalansky 
01CY  Attn  MNNXH  M  Delvecc';  ,.o 
01CY  Attn  MNN  W  Crabtree 
01CY  Attn  MNNXH  D  Gage 
01CY  Attn  MNNX 

Deputy  Chief  Staff 
Research,  Development  &  AC0 
Department  of  the  Air  Force 
Washington,  DC  20330 
01CY  Attn  AFRD0I  N  Alexandrow 
01CY  Attn  AFRDPN 
01CY  Attn  AFRDOI 

Commander 

Foreign  Technology  Division,  AFSC 
Wright-Patterson  AFB,  OH  45433 
0  ICY  Attn  NIIS  Library 

Director 

Strategic  Systems  Project  Office 
Department  of  the  Navy 
Washington,  DC  20376 
01CY  Attn  NSP-272 
01CY  Attn  NSP-43  (Tech  Lib) 

President 

Naval  War  College 

Newport ,  RI  02840 

01CY  Attn  Code  E-ll  (Tech  Service) 

Commanding  Officer 

Naval  Weapons  Evaluation  Facility 

Kirtland  Air  Force  Base 

Albuquerque,  NM  87117 

01CY  Attn  R  Hughes 

01CY  Attn  Code  10  (Tech  Lib) 
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Office  of  Naval  Research 

Arlington,  VA  22217 

01CY  Attn  Code  474  N  Perrone 


Superintendent 
Naval  Postgraduate  School 
Monterey,  CA  93940 
(Desires  No  CNWDI  Documents) 
01CY  Attn  Code  1424  Library 
01CY  Attn  G  Lindsay 


Commander 

Naval  Sea  Systems  Command 
Washington,  DC  20362 
01CY  Attn  SEA  09G53  (Lib) 
01CY  Attn  SEA  0351 


Commander 

David  Taylor  Naval  Ship  R&D  Ctr 
Bethesda,  MD  20084 
(CNWDI  only  Attn  Mrs,  M  Birkhead 
Code  5815.6) 

01CY  Attn  Code  L42-3  (Library) 


Commander 

Naval  Electronic  Systems  Command 
Washington,  DC  20360 
01CY  Attn  PME  117-21 

Headquarters 
Naval  Material  Command 
Washington,  DC  20360 
01CY  Attn  MAT  08T-22 


Commander 

U.S.  Army  Mobility  Equip  R&D  CMD 
Fort  Belvoir,  VA  22060 
(CNWDI  to  Army  MAT  Dev  &  Readiness  Com 
01CY  Attn  Drome-WC  (Tech  Lib) 


Commander 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 
01CY  Attn  Code  013  E  Cooper 
01CY  Attn  Code  4471  (Tech  Lib) 

Commanding  Officer 

Naval  Research  Laboratory 

Washington,  DC  20375 

(RD  &  RD/N  Attn  Code  1221  for  & 

FRD  Attn  Code  2628  FOR) 

01CY  Attn  Code  4040  J  Boris 
01CY  Attn  Code  2627  (Tech  Lib) 

01CY  Attn  4040  D  Book 

Officer  in  Charge 

Naval  Surface  Weapons  Center 

White  Oak  Laboratory 

Silver  Spring,  MD  20910 

01CY  Attn  P44  H  Galz 

01CY  Attn  Code  F31 

01CY  Attn  Code  X211  (Tech  Lib) 

Office  in  Charge 

Naval  Civil  Engineering  Laboratory 
Port  Hueneme,  CA  93041 
01CY  Attn  Code  L53  Forrest 
01CY  Attn  Code  L08A  Library 
01CY  Attn  Code  L51  J  Crawford 
01CY  Attn  L51  R  Murtha 

Commander 

Naval  Facilities  Engineering  Command 
Washington,  DC  20390 
01CY  Attn  Code  04B 

Commander 

U.S.  Army  Missile  Command 
Redstone  Arsenal,  AL  35898 
01CY  Attn  RSIC 
01CY  Attn  DRDMI-XS 

Commander 

U.S.  Army  Nuclear  &  Chemical  Agency 

7500  Baclclick  Road 

Building  2073 

Springfield,  VA  22150 

(Desires  only  1  cy  to  library) 

01CY  Attn  J  Simms 
01CY  Attn  Library 


Commandant 
US  Army  War  College 
Caplisle  Barracks,  PA  17013 
01CY  Attn  Library 

Commander 

US  Army  Foreign  Science  &  Tech  Ctr 
220  7th  Street,  NE 
Charolttesville,  VA  22901 
01CY  Attn  DRXST-SD 

Commander 

US  Army  Material  Dev  &  Readiness  CMD 

5001  Eisenhower  Avenue 

Alexandria,  VA  22333 

01CY  Attn  DRCDE-DL  Flynn 

01CY  Attn  DRXAM-TL  (Tech  Lib)uncl  only 

Director 

US  Army  construction  Energe  Res  Lab 
P.0.  Box  4005 
Champaign,  IL  61820 
01CY  Attn  Library 

Division  Engineer 

US  Army  Engineer  Div  Huntsville 

P.0.  Box  1600,  West  Station 

Huntsville,  AL  35807 

01CY  Attn  Hnded-SR 

01CY  Attn  Hnded-FD 

Commander 

Harry  Diamond  Laboratories 
Department  of  the  Army 
2800  Powder  Mill  Road 
Adelphi ,  MD  20783 

(CNWDI-Inner  Envelope-Attn:  Delhd- 
Rbh  FOR) 

01CY  Attn  DELHD-I-TL  (Tech  Lib) 

01CY  Attn  Chif  Div  20000 

Commander 

US  Array  Armament  Material  Readiness  Com 
Rock  Island,  IL  61202 
01CY  Attn  MA  Library 

Commander  and  Director 

US  Army  Cold  Region  Res  Engr  Lab 

P.0.  Box  282 

Hanover,  NH  03755 

01CY  Attn  Library 


Director 

US  Army  Eng  Waterways  Exper  Station 

P.O.  Box  631 

Vicksburg,  MS  39180 

01CY  Attn  J  Zelasko 

01 CY  Attn  Wessd  J  Jackson 

01CY  Attn  J  Strange 

01CY  Attn  Wesse  L  Ingram 

01CY  Attn  Library 

01CY  Attn  Wessa  W  Flathau 

Commander 

US  Army  Material  &  Mechanics  Rsch  Ct 
Watertown,  MA  02172 
(Address  CNWDI:  Attn  Document 
Control  FOR:) 

01CY  Attn  Technical  Library 
01CY  Attn  DRXMR-TE  R  Shea 
01CY  Attn  DRXMR  J  Mescall 

Commander 

US  Army  Concepts  Analysis  Agency 

8120  Woodmont  Avenue 

Bethesda,  MD  20014 

01CY  Attn  CSSA-ADL  (Tech  Lib) 


Commander 

US  Army  Engineering  Center 
Fort  Belvoir,  VA  22060 
01CY  Attn  Technical  Library 
01CY  Attn  ATZA 

Director 

US  Army  Ballistic  Research  Labs 
Aberdeen  Proving  Ground,  MD  21005 
01CY  Attn  DRDAR-BLV 
01CY  Attn  DRDAR-BLT  J  Keefer 
01CY  Attn  DRDAR-TSB-S  (Tech  Lib) 

Dep  Ch  of  Staff  FOR  OPS  &  Plans 
Department  of  the  Army 
Washington  DC  20310 
01CY  Attn  DAMO-NC 


31 


k> 


Director 

BMD  Advanced  Technology  Center 

Department  of  the  Army 

P.0.  Box  1500 

Huntsville,  AL  35807 

OICY  Attn  ATC-T 

01CY  Attn  ICRDABH-x 

01CY  Attn  ATC-T 


Chief  of  Engineers 

Department  of  the  Army 

Forrestal  Building 

Washington,  DC  20314 

01CY  Attn  DAEN-MCE-D 

0  ICY  Attn  DAEN-RDL 

01CY  Attn  DAEN-MPE-T  D  Reynolds 

Under  Secy  of  Def  For  Rsch  &  Engrg 
Department  of  Defense 
Washington,  DC  20301 
01CY  Attn  Strategic  &  Space  Sys 
Rm.  3E129 


Commander 
Field  Command 
Defense  Nuclear  Agency 
Kirtland  AFB,  NM  87115 
01CY  Attn  FCTMOF 
01CY  Attn  FCT 
01CY  Attn  FCPR 
01CY  Attn  FCTT 

Director  Joint  Strat  TGT  Planning 

0FFUTT  AFB 

Omaha,  NB  68113 

01CY  Attn  JLA 

01CY  Attn  DOXT 

01CY  Attn  XPES 

OICY  Attn  NRI-STINFO  Library 
01CY  ATtn  JLTW-2 


Commander 

BMD  Systems  Command 
Department  of  the  Army 
P.0.  Box  1500 
Huntsville,  AL  35807 
OICY  Attn  BMDSC-KW 
OICY  Attn  BMDSCH-HW  R  Dekalb 
OICY  Attn  BMDSCH  H  N  Hurst 
OICY  Attn  BMDSCH  RC  Webb 

Commandant 

NATO  School  (SHAPE) 

APO  New  York  09172 

OICY  Attn  US  Documents  Officer 


Chairman 

Department  of  Defense  Explo  Safety 
Board 

2461  Eisenhower  Avenue 
Alexandria,  VA  22331 
OICY  Attn  Chairman 

Chief 

Field  Command 
Defense  Nuclear  Agency 
Livermore  Branch 
P.0.  Box  808  L-317 
Livermore,  CA  94550 
OICY  Attn  FCPRL 


Horizons  Technology,  Inc. 

7830  Clairement  Mesa  Blvd 
San  Diego,  CA  92111 
OICY  Attn  R  Kruger 

Card,  Inc. 

7449  N  Natchez  Avenue 
Niles,  IL  60648 

OICY  Attn  G  Neidhardt  (unci  only) 


Staff 


I 
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Director 

US  Army  Materials  &  Mechanics 
Research  Center 
Watertown,  MA  02172 
01CY  Attn  J.  Dignam 

Commander 

US  Army  Air  Mobility  R&D 
Laboratory 

Fort  Eustis,  Va.  23604 
01CY  Attn  S.  Pociluyko 

Commander 

US  Army  Natick  R&D  Command 

Natick,  MA  01760 

01CY  Attn  DRDNA-U ST ,  J.  Siegel 

Director 

Waterways  Experiment  Station 
P.0.  Box  631 

Vicksburg,  Mississippi  39180 
0 ICY  Attn  W.  Flathau  (WESSE) 

Commander 

White  Sands  Missile  Range 

White  Sands  Missile  Range,  NM  88002 

01CY  Attn:  STEWS-TEAN,  J.  O'Kurna 

Commander 

US  Army  Missile  Command 

Rfedstone  Arsenal,  AL  35809 

01CY  Attn  DRSMI-RHB ,  Mr.  H.W.  Green 

Commander 

US  Army  Missile  Command 
Redstone  Arsenal,  AL  35808 
01CY  Attn  DRSMI-RRR ,  S.P.  Gibson 

Dept  of  the  Navy 

Theater  Nuclear  Warfare  Project 

Office  (PM  23) 

Washington,  D.  C.  20360 
01CY  Attn  Capt.  D.  Alderson 

Commander 

US  Chemical  Systems  Laboratory 
Aberdeen  Proving  Ground 
Aberdeen,  MD  21010 
01CY  Attn  DRDAR-CLW-E,  J.  Zuraw 


Director 

USA  Ballistic  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21019 

01CY  Attn  DRDAR-BLT,  D.  Menne 

01CY  Attn  W.J.  Taylor 

01CY  Attn  C.W.  Ktchens 

01CY  Attn  F.  Oertel 

0 ICY  Attn  F.  Oertel 

01CY  Attn  J.  Santiago 

01CY  Attn  J.H.  Keefer 

01 CY  Attn  N.  Ethridge 

01CY  Attn  W.  Schuman 

Director 

USA  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground,  Md.  21010 
01CY  DRDAR-BLV,  O.T.  Johnson 
01CY  R.  Vital 
01CY  D.  Rigotti 
01CY  D.  Haskell 

Assistant  Secretary  of  the  Navy 
(R  &  D) 

Department  of  the  Navy 
Washington,  D.C.  20350 

Chief,  Naval  Operations 
Department  of  the  Navy 
Washington,  D.C.  20350 
01CY  Attn  OP  985 

Chief 

Naval  Material  Office 
Department  of  the  Navy 
Washington,  D.C.  20350 
01CY  Attn  I.  Jaffee 


Commander 

Experimental  Engineer  Division 

Hanover,  NH  03755 

01CY  Attn  A.F.  Wuori  USA,  CRREL 

Director 

Nuclear  Program  Office 
Naval  Surface  Weapons  Center 
White  Oak,  MD  20910 
01CY  Attn  R.  Jenkins 
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Director 

Defense  Nuclear  Agency 
Washington,  D.  C.  20305 
01CY  Attn  Dr.  E.  Sevin 

Director 

Defense  Nuclear  Agency 
Washington,  D.C.  20305 
01CY  Attn  SPD  Col  M.  Kovel 
01CY  Attn  SPTD  T.  Kennedy 
01CY  Attn  SPTD  LTC  Flory 
01CY  Attn  SPAS  J.I.  Slaughter 

Commander,  Test  Command 
Defense  Nuclear  Agency 
Kirtland  AFB,  NM  87115 
01CY  Attn  FCT 

Federal  Emergency  Mangement 
Agency 

Washington,  D.C.  20472 
01CY  Attn  Dr.  D.A.  Bettge 

Defense  Intelligence  Agency 
Washington,  D.  C.  20301 
01CY  Attn:  DB4C2  (Carl  Wiehle) 

HODA 

Washington,  D.C.  20310 

01CY  Attn  DAMA/CSM/N  LTC  E.  Driscoll 

Dept,  of  the  Navy 
Theater  Nuclear  Warfare 
Project  Office  (PM  23) 

Washington,  DC  20360 
01CY  Attn  CAPT  D.  Alderson 


Commander 

US  Army  Materiel  Development  and 
Readiness  Command 
5001  Eisenhower  Avenue 
Alexandria,  Va.  22333 
01CY  Attn  DRCMD  B.  Dunetz 
0 ICY  Attn  DRCLDC  Dr.  H.  Gieske 
01CY  Attn  DRCDE-SSN,  D.  Bassett 

Office  of  the  Chief  of  Engineers 

Research  &  Development  Office 

Pulaski  Building 

20  Massachusetts  Ave,  N.W. 

Washington,  D.  C.  20314 

01CY  Attn  Dr.  James  Choromokos,  Jr. 

D0D  Explosives  Safety  Board 
Forrestal  Bldg.  6A145 
Washington,  D.  C.  20305 
01CY  Attn  T.A.  Zaker 

Director 

Defense  Advanced  Research  Project 
Agency 

1400  Wilson  Blvd. 

Arlington,  Va.  22209 

Commander 

USA  Health  Services  Command 
Fort  Sam  Houston,  Texas  78234 

Commander 

US  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709 
01CY  Attn  DRXRO-ZC 
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Commander 
USA  TRADOC 

Fort  Monroe,  Va.  21651 
01CY  Attn  ATCD-Z 

Commander 

USA  Corps  of  Engineers 
Huntsville  Division 
P.0.  Box  1600 
Huntsville,  A1  35807 

Commander 

US  Army  Nuclear  and  Chemical  Agency 
7500  Backlick  Road,  Bldg  2073 
Springfield,  Va  22150 
0 ICY  Attn  Col  L.E.  Minnich 
01CY  Attn  Maj  J.  Uecke 

Office  of  the  Chief  of  Engineers 
20  Massachusetts  Avenue,  N.W. 
Washington,  D.  C  20314 
0 ICY  Attn  DAEN/MPE-T  Mr.  D.S.  Reynolds 
01CY  Attn  DAEN/RDM  Mr.  J.  Healy 

Commander 

USA  Foreign  Science  and  Tech.  Center 
220  7th  Street,  N.E. 

Charlottesville,  Va  22901 
01CY  Attn  DRXST-SD  Charles  Ward 

Director 

Ballistic  Missile  Defense  Advanced 

Technology  Center 

P.0.  Box  1500 

Huntsville,  A1  35807 

01CY  Attn  M.  Whitfield 

01CY  Attn  M.  Capps  (ATC-T) 


Commander 

USA  Test  &  Evaluation  Command 
Aberdeen  Proving  Ground,  MD  21005 
01CY  Attn  DRSTE-SE  Mr.  R.  Galassso 

Commander 

USA  Armament  R&D 

Dover,  N.J.  07801 

01CY  Attn  DRDAR-LCW  Dr.  J.  Frasier 
Q1CY  Attn  DRDAR-LCN  P.  Angelotti 
0 ICY  Attn  DRDAR-LCS-V  M.A.  Ravotto 

Director 

AMSAA 

Aberdeen  Proving  Ground,  MD  21005 
0 ICY  Attn  DRXSY-S,  J.  Brand 

Commander 

US  Army  Electronics  Research  and 
Development  Command 
Fort  Monmouth,  N.J.  07703 
01CY  Attn  DRDEL-S  Dr.  W.  McAfee 

Commander 

CECOM 

Fort  Monmouth,  N.J.  07703 

01CY  Attn  DRSEL-C0M-RM-1,  A  Freer 

Project  Manager 
Pershing  Project  Office 
Redstone  Arsenal,  AL  35809 
•01CY  Attn  DRCPM-PE-EA  W.  Johnson 
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Commander 

Naval  Surface  Ueapons  Center 
White  Oak,  M d.  20910 
01CY  Attn  Code  241  M.M.  Swisdak,  Jr. 
0 ICY  Attn  Code  241  W.G.  Filler 

Director 

US  Army  Construction  Engineering 
Research  Laboratory 
P.0.  Box  4005 
Champaign,  II.  61820 
0  ICY  Attn  Dr.  W.E.  Fisher 
01CY  Attn  Alvin  Smith 


Commander 

Naval  Civil  Engineering  Laboratory 
Port  Hueneme,  Ca  93043 

Commander 

Naval  Weapons  Evaluation  Facility 
Kirtland  AFB,  NM  87117 
01CY  Attn  WEM( R.  Friedburg) 

Commander 

Naval  Sea  Systems 

Command 

Washington,  D.  C.  20350 
0 ICY  Attn  Code  SEA  322,  G.  Sullivan 

Director 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  Va.  22448 
01CY  Attn  J.  Yagla 

Commander 

Naval  Ship  Research  &  Development 
Center 

Bethesda,  Md.  20084 

01CT  Attn  Structures  Dept,  S.  Wang 


Naval  Research  Laboratory 
Lab  for  Computational  Physics 
Washington,  D.  C.  20375 
01CY  Attn  David  Book 
01CY  Attn  Mark  Fry 
01CY  Attn  David  Fyfe 
01CY  Attn  Pradeep  Kamat h 

Commander 

Naval  Coastal  System  Center 
Panama  City,  FI  32407 
01CY  Attn  R.W.  Buhl 

HO  US  Air  Force 

Washington,  D.  C.  20330 

01CY  Attn  RDOI,  Maj  N.  Alexandrov 

Commander 

Air  Force  Systems  Command 

Andrews  AFB,  MD  20331 

01CY  Attn  DLW ,  Maj.  R.  Padfield 

Commander 

Air  Force  Weapons  Laboratory 
Kirtland  AFB,  NM  87115 
0 ICY  Dr.  D.  Payton  (NT) 

01CY  M.A.  Plamondon  (DES) 

Commander 

Air  Force  Weapons  Laboratory 
Kirtland  AFB,  NM  87115 
01CY  Attn  R.  Guice  (NTED) 

01CY  Attn  DY 
01CY  Attn  A.  Sharp 
01CY  Attn  D.  Renick 

Commander 

USAF  Electronics  Systems  Division 
Hanscom  AFB,  Ma  01731 
01CY  Attn  R.  Karlson 

Commander 

Air  Force  Flight  Dynamics  Lab. 
Wright-Patterson  AFB,  OH  45433 
0 ICY  Attn  E.  Pugh 


Command 

Naval  Ship  Research  &  Development 
Center 

Portsmouth,  Va.  23709 

01CY  Attn  Underwater  Explosion  Research 
Division  Dr.  E.  Palmer 
01CY  Attn  R.  Higginbotham 


Sandia  National  Laboratory 
P.0.  Box  5800 

Albuquerque,  New  Mexico  87115 
01CY  Attn  M.  Vigil  (4522) 

0 ICY  Attn  J.  Reed  (1111) 

01CY  Attn  L.  Vortman 

Los  Alamos  Scientific  Laboratory 
P.O.  Box  1663 
Los  Alamos,  NM  87545 
0 ICY  Attn  Dr.  T.  Dowler 
0 ICY  Attn  Dr.  J.R.  .Lilley 

Lawrence  Livermore  Laboratory 

P.O.  Box  808 

Livermore,  Ca  94550 

01CY  Attn  H.  David  Glenn  (L-200) 

01CY  Attn  P.S.  Brown 

Oak  Ridge  National  Laboratories 
P.O.  Box  X 
Oak  Ridge,  TN  37830 
01CY  Attn  Health  Physics  Division 
C.K.  Kearny 


Aerospace  Corporation 
P.O.  Box  92957 
Los  Angeles,  Ca.  90009 
01CY  Attn  H.  Mirels 

Agbabian  Associates 
250  N.  Nash  Street 
El  Segundo,  Ca.  90245 
01CY  Attn  F.B.  Safford 

BDM  Corporation 
1920  Aline  Avenue 
Vienna,  Va.  22180 
01CY  Attn  W.  Sweeney 


HO  AFESC/RDCS 
Tyndall  AFB,  Flor 
0 ICY  Attn  Maj.  Ha 

USAF 

Armament  Laborato 
Eglin  AFB,  Florid 
01CY  Attn  P.T.  Na 

Commander 

Air  Force  Cambrid 

Laboratory 

L.G.  Hanscotn  Fiel 

Bedford,  MA  01730 

Commander 
US  Air  Force  Ball 
Corporation 
Norton  AFB,  Ca.  9 
0 ICY  Attn  MNNH  Ma 

Commander 
USAF  Aeronautical 
Wright-Patterson 
01CY  Attn  Materia 

US  Geological  Sur 
Geologic  Division 
2255  North  Gemini 
Flagstaff,  AZ  860 
01CY  Attn  Dr.  Dav 

Director  for  Mili 
US  Atomic  Energy 
Washington,  D.C. 

Sandia  National  L 
P.O.  Box  5800 
Albuquerque,  NM  8 
0 ICY  Attn  A.  Chab. 
01CY  Attn  W.R.  Da 
0 ICY  Attn  F.  Math' 


IIT  Research  Institute 
10  W  35th  Street 
Chicago,  II  60616 
01CY  Attn  D.  Anderson 
01CY  Attn  Keith  McKee 

Jaycor 

205  S.  Whiting  Street 
Alexandria,  VA  22304 

Kaman  Sciences  Corp. 

1911  Jefferson  Davis  Highway 
Arlington,  VA  22202 
01CY  Att  Dr.  D.C.  Sachs 

Kaman  Sciences  Corporation 
1500  Carden  of  Cods  Road 
P.0.  Box  7463 

Colorado  Springs,  CO  80933 
01CY  Attn  G.L.  Roark 

Kaman  Tempo 
7800  Marble  Ave. ,  NE 
Albuquerque,  NM  87110 
01CY  Attn  Robert  H.  Miller 

Kaman  Tempo 
715  Shamrock  Rd. 

BelAir ,  MD  21014 
01CY  Attn  E.J.  Bryant 

Kaman  Tempo 

Huntington  Bldg.  Suit  506 
2560  Huntington  Ave 
Alexandria,  VA  22303 
01CY  Attn  J.  Petes 


Boeing  Aerospace  Co. 

P.0.  Box  3999 
Seattle,  WA  98124 

01CY  Attn  R.M.  Schmidt  (M/S  42-37) 

Calspan  Corp/ATC 
P.0.  Box  400 
Buffalo,  NY  14225 
01CY  Attn  M.G.  Dunn 

Civil  Systems  Incorporated 
P.0.  Box  2083 
Midland,  TX  79702 
01CY  Attn  Mr.  Steve  Melzer 

Consolidated  Developments,  Inc. 

515  Sherwood  Drive 
Marion,  VA  24354 
01CY  Attn  David  E.  Cary 

Electro-Mechanical  Systems  of 
Albuquerque 
P.0.  Box  11730 
Albuquerque,  NM  87192 
01CY  Attn  Dr.  Ray  A.  Shunk 

H-TECH  Laboratories,  Inc. 

P.0.  Box  1686 

Santa  Monica,  CA  90406 

01CY  Attn  B.A.  Hartenbaum 

Information  Science,  Inc. 

123  West  Padre  Street 

Santa  Barbara,  CA  93105 

01CY  Attn  Dr.  Walter  F.  Dudziak 
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Dr.  Marion  Balcerzak 
GATX 

7449  Matchez  Ave 
Niles,  ILL  60648 

Mr.  Burton  S.  Chambers 
SEA,  Inc. 

P.O.  Box  308 
McLean,  VA  22101 

Burt  H.  Collins 
Tech.  Reps.  Inc. 

5000  Marble,  NE,  Suite  222 
Albuquerque,  NM  87110 

Dr.  Steve  Gill 
ARTEC  Associates 
26046  Eden  Landing  Road 
Hayward,  CA  94545 

Dr.  J.P.  McLain 

New  Mexico  Inst,  of  Mining  &  Technology 
Socorro,  NM  87801 

J.  Kelso 

7124  Wolf tree  Lane 
Rockville,  MD  20852 


SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 
01CY  Attn  Mr.  S.B.  Martin 
01CY  Attn  J.D.  Colton 

Southwest  Research  Institute 
P.O.  Box  28510 
San  Antonio,  TX  78284 
01CY  Dr.  W.  Baker 

Scientific  Service,  Inc. 

517  E.  Bayshore  Road 
Redwood  City,  CA  94063 
01CY  Attn  Dr.  B.  Gabrielsen 
01CY  Attn  C.  Wilton 

Systems,  Science  and  Software 
P.O.  Box  8243 
Albuquerque,  NM  87198 
01CY  Attn  C.E.  Needham 

Systems,  Science  and  Software 
P.O.  Box  1620 
La  Jolla,  CA  92038 
01CY  Attn  T.H.  Pierce 
01CY  Attn  R.E.  Duff 


Higgins,  Auld  &  Associates,  Inc. 
2601  Wyoming  Blvd.  NE  Suite  H-l 
Albuquerque,  NM  87112 
0ICY  Attn  Dr.  C.J.  Higgins 
01CY  Attn  J.  Bratton 

Dr.  Donald  R.  Richmond 
Lovelace  Biomedical  &  Environmental 
Research  Institute 
P.O.  Box  5890 
Albuquerque,  NM  87115 


Temet  USA,  Inc. 

P.O.  Box  761 

Vienna,  VA  22180 

01CY  Attn  Dale  K.  Henderson 

TWR  Electronics  &  Defense 
1  Space  Park 
Redondo  Beach,  CA  90278 
01CY  Attn  D.R.  Ausherman 
01CY  Attn  J.J.  Schrader 


Kendan  Associates 

801  E.  Charleston  Rd.,  M 

Palo  Alto,  CA  94303 

01CY  Attn  Robert  E.  Dannenbert 


Lockheed  Missile  and  Space  Co. 
P.0.  Box  504 
Sunnyvale,  CA  94088- 
01CY  Attn  M.  C.  Vivian 

Los  Alamos  Technical  Assoc,  Inc. 

P.0.  Box  410 

1650  Trinity  Drive 

Los  Alamos,  NM  87544 

01CY  Attn  P.  Hughes 


Management  Science  Assoc. 

P.0.  Box  239 

Los  Altos,  CA  94022 

01CY  Attn  K.  Kaplan 

Marco  Scientific 
1031-H  East  Duane  Ave 
Sunnyvale,  CA  94086 
01CY  Attn  J.H.  Owren 

MITRE  Corporation 

P.0.  Box  208 

Bedford,  Mass  01730 

01CY  Attn  Mr.  Edward  M.  Sullivan 

McDonnel  Douglas  Astronautics  Co. 
5301  Bolsa  Ave. 

Huntington  Beach,  CA  92647 
01CY  Attn  Dr.  P.  Lewis 

Martin  Marietta  Corporation 
P.0.  Box  5837 
Orlando,  FL  32855 
01CY  Attn  G.  Fotieo 


Director 

New  Mexico  Engineering 
Research  Institute 

University  of  New  Mexico 
Albuquerque,  NM  87131 
01CY  Attn  G.  Leigh 
01CY  Attn  H.W.  Wampler 


Physics  Applications,  Inc. 
2340  Harris  Way 
San  Jose,  CA  95131 
01CY  Attn  Coye  Vincent 

Physics  International 
2700  Merced  Street 
San  Leandro,  CA  94577 
01CY  Attn  F.M.  Sayer 
01CY  Attn  J.M.  Thomsen 


Pacific-Sierra  Research  Corp. 
1456  Cloverfield  Blvd. 

Santa  Monica,  CA  90404 
01CY  Attn  H.  Brode 

Rand  Corporation 
1700  Main  Street 
Santa  Monica,  CA  90406 
01CY  Attn  Dr.  0.  Nance 

San  Diego  State  University 
Dept,  of  Civil  Engineering,  SDSU 
San  Diego,  CA  92182 
01CY  Attn  Dr.  R.J.  Fragasgy 

Martin  Marietta  Aerospace 
P.0.  Box  179 
Denver,  COLO  80201 
01CY  Attn  Robert  G.  DeRaad 


